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Abstract 



Starting from the SO{2, 2n) Chern-Simons form in (2n + 1) dimen- 
sions we calculate the variation of conserved quantities in Lovelock gravity 
^ ' and Lovelock-Maxwell gravity through the covariant formalism developed 

0^ , in [23]. Despite the technical complexity of the Lovelock Lagrangian we 

obtain a remarkably simple expression for the variation of the charges en- 
suing from the diffeomorphism covariance of the theory. The viability of 
the result is tested in specific applications and the formal expression for 



o 

00 

' the entropy of Lovelock black holes is recovered. 
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5-1 '■ 1 Introduction 

Lovelock gravity can be considered as the most general extension of General 
Relativity to higher dimensions: in fact, even if it contains higher curvature in- 
teraction terms, it is nevertheless governed by field equations which are just of 
' the second order in the metric field [39, 40]. Far from being just a mathematical 

curiosity. Lovelock gravity as well as higher curvature theories in general, have 
received in the recent past a renewed interest motivated by the hope of learning 
something about the nature of quantum gravity. As a matter of facts, higher 
curvature interaction terms arise in the study of back-reactions of quantum field 
energy as well as in the low-energy limit of string theory (see [9, 11, 30, 42]). 

In this paper we address the problem to calculate the charges associated 
to the diffeomorphism symmetries of Lovelock gravity in any odd dimensional 
spacetime. To this end we make use of the recipe worked out in [23] to al- 
gorithmically define (the variation of) conserved quantities directly from the 
equations of motion. One of the major advantages of the formalism developed 

* allemandi@dm.unito.it 
tfrancaviglia@dm.unito.it 
traiteri@dm.unito.it 



1 



in [23] is that it enables to establish the transformation rules acting on conserved 
quantities in the transition from different, but nevertheless on-shell equivalent, 
Lagrangian theories. This property is well suited to tackle with our problem. It 
was indeed shown in [6, 13] that a special case of Lovelock gravity (with a fixed 
choice of the Lagrangian coefficients) in D = 2n + 1 dimensions can be derived 
from a Chern-Simon theory of the group SO{2,2n) (or SO{l,2n + 1)). Com- 
bining together this remarkable result with the formalism of [23] for conserved 
quantities, we shall show that the calculation of conserved quantities for Love- 
lock gravity becomes rather easier if we first compute them in the Chern-Simons 
framework and we take then into account the suitable transformation rules on 
the field variables as well as on the infinitesimal generators of symmetries. Even 
though we end up in this way with an expression which holds true only for a very 
special case of Lovelock gravity, the contribution to conserved quantities coming 
from each higher curvature interaction term can be nevertheless isolated. This 
property allows one to derive a general formula for (the variation of) conserved 
quantities for any Lovelock theory, i.e. with completely arbitrary choice of the 
coefficients in front of each curvature interaction term. 



Let us now analyse the problem in detail. The D = 2n+l dimensional Chern- 
Simons Lagrangian for the anti-de Sitter gauge group 50(2, 2n) is written in 
terms of a gauge connection 

A = A^^dx^JAB (1) 

where = 0, . . . , 2n; A, B = 0, . . . , 2n + 1; and Jab are the Lie algebra genera- 
tors: 

[Jab, Jcd] = {-JacVbd + JadVbc + JbcVad - JbdVac} (2) 

with r]AB = (-1, 1, • • • , 1, -1). 

The Lagrangian turns out to be [13]: ^ 

(ji A) =K J^dt<AA {tdA + eA^Y > (3) 
where < > denotes the invariant form: 

< JAiBi ■ ■ ■ JAn+lBn+l >= ^AiBi...A„ + iBn + l (4) 

According to [6] the constant k in (3) is fixed as: 

^= (D-2)!Ob-2 ' ^ = 

where 0^-2 is the area of the (D — 2) dimensional unit sphere and I has the 
dimension of a length. The Euler-Lagrange field equations ensuing from (3) 
turn out to be: 

^AAit^B^^. - ^eA,B ,...A„^,B.^,F^-^^ A ... A i^^"^" = (6) 

^In the sequel we shall make use of the notation f{j^y) to denote a quantity " / " depending 
on the variables y together with their derivatives up to the order k. 
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where F^^ = dA"^^ + A^c f\ A^^ is the field strength of the gauge connec- 
tion. Furthermore, the set of generators Jab can be split as Jab = {JabyPa), 
a = 0, — 1; namely: into the generators Jab of Lorentz rotations in D di- 

mensions and the generators Pa := Ja,2n+i of inner translations (see for example 
[1], [13], [51]). Accordingly, the SO{2,2n) gauge connection can be decomposed 
as follows: 

A = Lo'^'jab + ye«P„ (7) 
This assumption implies that the field strength in turn decomposes as: 



where 



F = (i?"'' + -2 e'' A e") J„b + -T'^ Pa (8) 



j^ab ^ ^ab ^ ^a^ ^ ^cb (9) 

is the field strength of the •50(1, 2n) gauge connection w"'', while 

T" = de" + uj^b A (10) 

is the torsion 2-form. 

Substituting expressions (7) and (8) into (3) and taking into account that 
the only terms which do not vanish are those of the kind < Ja^bi ■ ■ ■ Ja^b^Pa >= 
eai6i...a„b„a we Can rewrite the Chern-Simons Lagrangian in terms of the new 
dynamical fields (w,e). Apart from boundary terms we obtain: 

n 

L2{e,fLo) = K^CpCp{e,fu;) (11) 

where 

"""^ {D- 2p) {p) ^^^^ 

Cp{e,j^Lj) = ea^...aoR"^"^ A ... A i?«2p-ia2p ^ ga2p+i A ... A 6°^ (13) 

The Lagrangian (11) turns out to be the smu of dimensionally continued Euler 
densities [6, 13, 47] and we shall refer to it as the Euler-Chern-Simons La- 
grangian. Notice that in the expressions (13) the "vielbein" e" appears without 
derivatives and the gauge connection o;"^ enters only through its field strength 
R. This remarkable property is due to the fact that all the terms involving 
explicitly the torsion and/or the gauge connection w are pushed, through inte- 
grations by parts, into boundary terms which are discarded in the expression 
(11). Notice also that, until field equations are not yet solved, no a priori rule 
exists relating the fields e" and w"*": they are completely independent.^ 

Geometrically speaking we are assuming that the Lagrangian L2 is based on a con- 
figuration bundle Y which is a product bundle Y = Y\ Xm Y2 with fibered coordinates 
(a;'',eS,a;^=), where {xi^,el) and (x'^ , Wp^) are fibered coordinates on Y\ and Y2, respectively. 
The configuration bundle Y is thence a gauge natural bundle admitting a principal SO(l, 2n) 
bundle Q as structure bundle [18, 38]. Accordingly Yi is the spin-frame bundle [19] while 
Y2 = J^Q/SO{l,2n) is the connection bundle; see [38, 45] 
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Varying the Lagrangian (11) with respect to the independent fields {e,w) we 
obtain, respectively, the field equations: 



Se" — I '^ooiOi-'-anOn-'^ /\ . . . /\ n — u 

= nJlfz . 1, T<* A A A pan-ibn-i = n 

i5ti, <>„()„ — ; taaiOi...a„o,i-' /\ -fi, /\ . . . /\ — u 



(14) 



where 



Notice that field equations (14) are dynamically equivalent to field equations (6) 
once the substitution (7) is taken into account. Namely, the difference between 
the two sets of equations is just a matter of notation. 

If we now consider the particular solution = of the second set of field 
equations (14), it turns out that the 50(1,272) connection u)"^ is the spin con- 
nection of the vielbein e°: 

a;V = e«(S«^^-I]^% + S^;3«)ef, SV = e?5[^e^] (16) 

and 

R'^\. = e^e'^RP^^ifg) (17) 

where R'^'^ g) is the Riemann tensor of the metric = rjabelj^et- 
serting (17) back into the Lagrangian (11) we end up with a special Lovelock 
Lagrangian, namely: 

n 

Laifg) := L2{e,j^oj)\^^^^ = KY,CpCp{fg) (18) 

p=0 



where: 



\ (2n-2p)! 

J lD-2p 



(19) 



^Af9) = ^V9 siii:::l^:R^2 ■ ■ ■ ^ZV-Zl (20) 

and Svlllll^^p denotes the totally skew-symmetric product of Kronecker deltas 
normalized to have values 0, ±1.*^ In particular: Cq = gives rise to the term 
corresponding to a cosmological constant A proportional to (or propor- 

tional to 1/P if we had started with the group SO{l, 2n+ 1)); £i = y/g R is the 
Hilbert Lagrangian in dimension 2n+l and £2 = (-^^~'^-^|-^a+-^a/3-^/iJ^) 
the Gauss-Bonnet term. Notice that the coefficients in the Lovelock Lagrangian 
(18) are uniquely fixed and in fact dictated by the initial Chern-Simons form 

^In writing (20) we made use of the fact that, according to our notation, 

<:M1 •■•M2pC«2p + l- -"D <: — ( n — ^nV ri'^^'"'^^'' (9^\ 



4 



(3); sec [6, 13, 15]. From now on we shall refer to the Lagrangian (18) as 
Lovelock-C'hern-Simons Lagrangian. 

Let us now siimmarize. We have introduced so far three different La- 
grangians; see expressions (3), (11) and (18). The transition among them can 
be schematically drawn as follows: 

L.ifA) ^-^-'^ L2(e,/a;) ^ L,{fg) (22) 

Notice however that the equivalence between the first two Lagrangians, as al- 
ready explained, holds true off-shell. On the contrary, the equivalence between 
Z/2 and L3 is just an on- shell equivalence since it holds true only along the 
space of solutions with T — Q. 

The goal of the present paper is to calculate the conserved quantities, in 
particular the energy, for each one of the Lagrangians in (22). We point out 
that this task is more subtle than it could first look like. The reason is mainly 
due to two different causes. 

1. Let us suppose that we are somehow able to calculate the charge Q{Li, E) 
associated to the Lagrangian Li{j^A) and relative to an infinitesimal gen- 
erator of symmetries S which is common to both the Lagrangians Li 
and L2. According to the prescription (22) we could be induced to iden- 
tify the charge Q{L2,'S), relative to the second Lagrangian in (22) and 
associated with the same generator S, with the expression Q{L2,'E) = 
Q{Li,'E.)\^^^ i.e. by merely inserting the splitting (7) into (5(Li,S). 
This rule could fail in giving the right expression. Indeed we stress that 
in the transition from the Lagrangian (3) to the Lagrangian (11) diver- 
gence terms have to be discarded. As far as field equations are concerned, 
those divergence terms are clearly irrelevant. On the contrary, they can 
become strongly relevant when dealing with conserved quantities. For ex- 
ample, conserved quantities calculated via Noether theorem, Hamiltonian 
or Hamilton-Jacobi-like methods are all sensitive to divergence terms of 
the Lagrangian (see, e. g., [10, 12, 24, 32, 33, 44]). 

2. Let us now suppose we have been somehow able to bypass the afore- 
mentioned problem. There remain still another open question. Indeed, 
conserved quantities of a field theory are mathematical expressions which 
are attributed the meaning of physical observables. They are built out by 
keeping into account the symmetry properties of the field theory we are 
dealing with. Hence, we have to know the relationships interplaying among 
the symmetries of the three Lagrangians (22) (or their field equations) if 
we want to establish a correspondence among their charges, too. Entering 
into the details of the matter, we observe that the Chern-Simons theory 
is a gauge natural theory; see [18, 19, 20, 27, 38]. It is based on a config- 
uration bundle C which is the bundle of connections C = P/ SO{2,2n), 
where P denotes a principal bundle with Lie group SO{2,2n), called the 
structure bundle of the theory; see [45]. Fibered coordinates on C are 
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The most general infinitesimal generator of symmetries for 
the equations of motion (6) is a projectable vector field Sp on P of the 
kind 

Sp=e5^ + S^^PAB (23) 

where pab is a local basis of right invariant vector fields on P (in a triv- 
ialization {x,g) of P we have p = (gd/dg), g G SO{2,2n)). The vector 
field Sp induces a transformation of the dynamical fields according to the 
rule: 

£e,A^'' = i'd.Afi^ + d^eAf^+ D, (24) 

and pure gauge transformations correspond to ^'^ = 0. Therefore the 
charge Q{Li,Ep) associated to this symmetry must depend, in general, on 
the dynamical field A, on the coefficients and together with their 
derivatives up to a fixed (finite) order. If we now consider the splitting 

S^^ = (S°^S°) (25) 

of the vertical part of the vector field Sp into "rotational" part and "trans- 
lational" part, we obtain from (7) and (24) 

£^pel = ^Pd,e; + d^ee; + l D^E- -E^e"^ (26) 

and, if we do not assume a priori the condition = 0, field equations 
(14) are separately invariant for the transformation (26) only when we set 
E" = 0> 

It follows that the most general infinitesimal generator of symmetries for 
each one of the equations of motion (14) is a projectable vector field Eq 
on the relevant SO{2n, 1) principal bundle Q of the theory (see footnote 
2) which, in turn, is a subbundle of P. In coordinates: 

SQ = e''5^ + S«Vb (27) 

where now pat is a local basis of right invariant vector fields on Q. The 
charge Q{L2, Eq) associated to the vector field (27) must depend, in gen- 
eral, on the dynamical fields (e, lo), on the coefficients and S**^ together 
with their derivatives up to a fixed order and can be suitable obtained 
from Q{Li,Ep) through the prescription (7) on the dynamical variables 
and the prescription = on the generators of symmetries. 



^Basically, the Chern-Simons symmetry is broken when we declare that the dynamical 
fields {e,Lj) are independent fields or, in other words, when we assume that the configuration 
bundle Y = Yi Xm ^2 has a product bundle structure. From (26) it is clear that "inner 
translations" of the AdS group, namely those generated by the components S", induce trans- 
formations on the configuration bundle which do not respect the product bundle structure. 
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Let us finally draw our attention to the Lovelock Clicrn- Simons Lagrangian 
(18). The latter Lagrangian is a natural (also said covariant or diffeomor- 
phism invariant; see [25, 33, 38]) Lagrangian and the generators of sym- 
metries are (not necessarily Killing) vector fields ^ = ^^9^^ on the base 
manifold (i.e. the D-dimensional spacetime). It follows that the charges 
Q{L3,(^) must depend, in general, on the dynamical field g and on the 
coefficients together with their derivatives up to a fixed order. To ob- 
tain a correspondence between Q{L2,'B.q) and Q{L3,^) the prescriptions 
(16) and (17) alone are not sufficient since we also need a prescription 
to fix the coefficients S"'' = $'^''(^^, rfiy^^, • • •) in terms of the coefficients 
^'^ and their derivatives up to a fixed order, namely we need to select a 
lift of the vector field ^ up to the structure bundle Q of the Lagrangian 
L2. Notice that the rule S"'' = is only local and it is not preserved (in 
general) by the automorphisms of the structure bundle Q, meaning that 
we have to resort to some other covariant lift. If we are able to select a 
lift S"'' = $''''(^^, (iy^^, . . .) which is mathematically well-defined (i.e. a 
covariant lift) as well as physically viable (i.e. it allows to reproduce the 
expected physical values for conserved quantities in specific applications) 
we have the following scheme: 

Q{L^{j^A),Ep) ^=p>-^^ Q{L2{e,j'u,),EQ) ^ Q{Ls{g),^) 

(28) 

governing the transition among conserved quantities. 



In the sequel we shall follow a recipe to calculate conserved quantities which 
has been recently worked out in [23] . The formalism developed in [23] is in fact 
well suited to tackle the issue we are dealing with, mainly because it allows to 
overcome both the problems we have just outlined. Basically, the construction 
of conserved quantities is defined as follows. Given a field theory the dynamics 
of which is derived from a Lagrangian L and given a vector field S which is an 
infinitesimal generator of symmetries for the equations of motions ^ , one can 
construct the variational Lagrangian L'{L,'E) through the contraction of field 
equations with the Lie derivatives of fields, i.e. L'{L, S) = — ^i^Hj/'. Handling 
the variational Lagrangian with the tools of Variational Calculus in jet bundles 
it is possible to algorithmically define the variation 5Q{L' ,'E) of covariantly 
conserved quantities along a one-parameter family of solutions. 

Since the variation of conserved quantities is obtained directly from the 
equations of motion the drawback related to problem 1) above is automatically 
overcome: the procedure is not affected by whichever divergence term we add to 
a given Lagrangian since field equations are insensitive to it. Moreover, also the 
problem raised in problem 2) can be easily bypassed. Indeed, by imposing the 
equivalence between the variational Lagrangian _L'(L2,2q) (built out starting 
from the equations of motion (14) and the generator of symmetries (27)) and 
the variational Lagrangian L'(L3,^) (where ^ is any spacetime vector field) a 
preferred lift Sq = $(^'',dy^'', . . .) of spacetime vector fields is automatically 
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ruled out from the theory and it allows to complete the diagram (28). The se- 
lected lift turns out to be the so-called generalised Kosmann lift; see [19, 27, 28]. 

Once we have imdcrstood the rules which allow us to make the transition 
from the (variation of the) conserved quantities 6Q{L[,'E.p) ensuing from the 
Chern-Simons Lagrangian (3) and the (variation of the) conserved quantities 
SQ{L'^,£^) associated to the Lagrangian (18), the latter quantities can be easily 
computed from the former ones. Moreover, even if the quantities obtained in 
this way are relative to the whole Lagrangian (18), the contribution coming 
from each Lagrangian Cp can be easily recognized and, therefore, we can define 
the conserved quantities for a generic Lovelock Lagrangian with arbitrary co- 
efficients in front of each term. Hence, despite the technical complexity of the 
Lovelock Lagrangian, we obtain a remarkably simple expression for its charges. 



The rest of the paper is organized as follows. In section 2 we shall shortly 

summarize the formalism developed in [23] to obtain the variation of conserved 
quantities directly from the equations of motion. We then apply it to calcu- 
late conserved quantities relative to the Chern-Simons form (3). Following the 
recipe schematically drawn in (28) we are then able to define the variation of 
conserved quantities for Lovelock gravity. A straightforward generalization of 
the aforementioned formalism will allow the generalization to Lovelock-Maxwell 
gravity. In section 3 the viability of the results previously obtained is tested for 
a five-dimensional charged black hole solution and the first law of black hole 
mechanics directly follows from the homological properties conserved quantities 
must satisfy. In section 4 the formal expression for the entropy of Lovelock sta- 
tionary black holes is finally calculated and it agrees with previous definitions; 
see [33, 35, 49]. 

2 Conserved Quantities 

In order to make the paper self-contained we shall here summarize the formalism 

to obtain the variation of conserved quantities directly from the equations of 
motion. Since we are interested only in the application of the technique to 
the Chern-Simons theory and to Lovelock gravity, we shall skip the rigorous 
geometric framework and we refer the reader to [23] for a deeper insight into 
the mathematical details of the matter. To stimulate the interest of physically- 
oriented readers we shall keep as much as possible the formal setting and the 
technical details to a minimum. 

Let us consider a field theory described through a Lagrangian of order k. It 
can be locally written as L = £(y', . . . j^^^..^^) where £ is the Lagrangian 
density, ds = dx^ A ... A dxP is the local volume element of the D-dimensional 
spacetime and we have collectively denoted with all the dynamical fields. 

The variation dxjO. of the Lagrangian density with respect to a vector field 
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X = = Sy^ can be generally written, through a well known integration 

by parts procedure, as follows: 

SxCij^'y) = ei{j^''y)X' + dxF^{f^-^y,j^-^X) (29) 

or, in terms of differential forms, as 

5xL = e{L, X) + dF{L, X) (30) 

Here the D-form 

e{L,X) = ei{f''y)X'ds (31) 

and the {D — l)-form 

F{L,X) = F\f''-^y,f-'X)dsx, dsx = h\ds (32) 

are called, respectively, the Euler-Lagrange form and the Poincare-Cartan form. 
Obviously field solutions are those field configurations = cr^{x) which satisfy 
the Euler-Lagrange equations: 

e.(/'=y)|,_(,)=0 (33) 

Let us now suppose that S is an infinitesimal generator of symmetries, namely a 
projectable vector field S = ^^{x)-^ + S'(a;, y)-^ which leaves field equations 
invariant: 

feei(j2'=j/)=0 (34) 

The contraction of the field equations (33) with the Lie derivatives £Ey^ of the 
field originates, apart for a sign, the so-called variational Lagrangian 

L'{L,E):=-ei{f''y)£sy'ds (35) 

which depends on the dynamical fields and the components and of 
the vector field S; it is clearly vanishing along solutions. According to the first 
variation formula (30), the variation of the variational Lagrangian L' yields:^ 

6xL' = e{L', X) + dF{L', X) (36) 

Moreover, since we have assumed that S is an infinitesimal generator of sym- 
metries it turns out that the variational Lagrangian is a pure divergence and 
therefore it is variationally trivial; see [23]. This means that e{L',X) = and, 
consequently, (36) reduces to the conservation law 

dF{L',X) = SxL' c^O (37) 

where ~ denotes equality on-shell. Moreover the Poincare-Cartan morphism 
F{L', X) ill (37) is linear in the Lie derivatives £sy and their derivatives up to 

^We remark that, while it is not restrictive to set S^i^ = as the components do not 
depend on the dynamical fields, we shall instead allow ^ 0. 
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a fixed order. Since in all physically relevant theories (e.g. in all gauge natural 
theories; see [18, 20, 38]) the Lie derivatives of the fields with respect to a vector 
field can be written as linear combinations in the components of the vector 
field itself together with their (covariant) derivatives, it turns out that also the 
Poincare-Cartan morphism F{L', X) is a linear combination of the coefficients 
and of the vector field S together with their derivatives up to a finite 
order which depends on the particular theory we are dealing with. This allows 
us to implement the so-called Spencer cohomology [18, 20, 25, 29]. Namely, the 
Poincare-Cartan morphism F{L',X), through repeated covariant integrations 
by parts with respect to the components of S, can be alghoritmically rewritten 
as [23]: 

F{L', X) = F{L', X) + dU{L', X) ~ dU{L', X) (38) 

where the term F{L',X) turns out to be proportional to the field equations 

together with their derivatives and is thence vanishing on shell, while the (£) — 2) 
form U{L' , X) is the (covariant) potential. The potential depends in turn on the 
dynamical fields y*, on the components (^'*, S*) of the infinitesimal generator of 
symmetries and on the components of the variational field X together with their 
derivatives up to a suitable order. Given a field configuration y = (t{x) and a 
{D — l)-dimensional region S of spacetime with boundary dT,, the integral 

6xQ^{<J,^) := I F{L',X) = / F{L',X) + [ U{L',X) (39) 

~ / U{L',X) (40) 

is assumed to be the variation 5xQ of the charge Q along a one-parameter 
family of field configuration described by X and relative to the symmetry origi- 
nated by !B. If y = a{x) is a solution of field equations and X is a solution of the 
linearized field equations (i.e. X is tangent to the space of solutions) the quan- 
tity SxQn (cr, S) ~ Jgj, U {L\ X) describes the variation of the charge Q moving 
from a to nearby solutions along a path in the space of solutions determined 
by X. Moreover if the region S is a (portion of) a Cauchy hypersurface with a 
boundary 9E and, in addition, the spacetime projection ^ of the vector field S 
is transverse to S, we agree to identify SxQt. with the variation of the energy 
E{dT:,l;, X) enclosed inside 3E. We stress that, once we fix 9S, there exist 
many ■'energies", depending both on the symmetry vector field ^ and on the 
variational vector field X; [2, 12, 26, 37]. The former vector determines, via its 
fiow parameter, how the hypersurface S evolves as the "time" flows. The vec- 
tor field X fixes instead the boundary conditions on the dynamical fields (e.g. 
Dirichlet or Newmann boundary conditions). Each choice of the pair {^,X) 
gives rise to a particular realization of a physical system endowed thence with 
its own "energy content" . 

Summarizing, we obtained the variation of conserved quantities starting from 
the equations of motion through a two-step procedure. The first step is nothing 
but an integration by parts of the variational Lagrangian with respect to the 
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vector field X (sec (36)). The second step amounts to perform a second inte- 
gration by parts of the Poincare-Cartan morphism F{L', X) with respect to the 
vector field S (see (38)). We remark that, a priori, there exists no guarantee 
that the resulting quantity SxQi: is in fact the variation of a functional Q of the 
intrinsic geometry of the boundary dT,. This becomes true only if appropriate 
boundary conditions are chosen; see [4, 8, 10, 12, 14, 26, 32, 37, 44]. 

The viability of such a procedure has been tested in [23] in the realm of 
natural and gauge natural theories of gravitations. The physically expected val- 
ues (or better, the commonly accepted expressions) for conserved quantities, in 
particular for the energy, have been there reproduced.^ Encouraged by those 
positive results we aim now to apply the formalism to Chern-Simons and Love- 
lock theories. 



The variational Lagrangian L'{Li, S) for the Chern-Simons theory is easily 
obtained from the equations of motion (6) : 

L'{Li,E) = -KeA,B,...A„^A^,F^'''' A ... A f^"B„ ^^^^„+iB„+i ^^^^ 

Taking into account that the Lie derivative (24) can be conveniently rewritten 
in a manifestly covariant form as: 

(A) 

£^Ai^ = ^Ff,^ + D X Sf,f ( Sf^f = + A^^n (42) 
the variation of (41) turns out to be (first step): 

5xL'{L^,E) = d,jr«(L'(Li,H),X) (43) 

where 

^"(L'(ii,S),X) = -|^e^^'^'-'^-'--'^"-6ABCi.A,B,...A„_,s„_, 

X {2n5'^^5At^£^A'i^+F^^^ (2€[«<5^1<5^^^ + <5,«<5Sf^f ) } 

Inserting (42) into (44) and integrating by parts the term containing the covari- 
ant derivative we obtain (second step): 

J^«(L'(Li,S),X) =;r«(L'(Li,S),X) + d/3ZY«^(L'(Li,S),X) (45) 

where 

^"(i'(Ll,S),X)=<5{-|^e"''--••'^"-eABA,B,...A„B„i^,t.f^...F/J.^ 

(46) 

^In the application [23, 26] of the formalism to staxdaxd gravity we found that the Brown- 
York quasilocal energy [12] can be obtained when ^ = 9( is the vector generating an ADM 
foliation of spacetime and X corresponds to Dirichlet boundary conditions of the induced 
boundary three-metric. 
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vanishes on-shell, while 

are the coefficients of the covariant potential (2n — l)-form U = ^U^^dSaff 
with dsais = ip\ ia\ ds (an expression similar to (47) was also obtained in [46] 
through a diffcrc;nt technique;). Wo remark that the integration of the potential 
gives rise, through expression (40), to the variation of the charge associated with 
the infinitesimal symmetry S. Notice that the potential can be rewritten in a 
more compact form as follows^: 

l^{L'{Li,E),X) = UK eABCDAiBi...An-iBn-i 

xF^i^i A ... A F^"-iBn-i ^ M^^S^^ (49) 

Since the potential has been obtained directly from the equations of motion 

and the equations of motion (6) are equivalent to field equations (14) we can im- 
mediately obtain the potential U{L'{L2,S),X) relative to the latter equations. 
According to the schema (28) it is given by: 

U{L' {L2, X) — —€abcaibi...a„-ibn-i 

A ... A R'^-^bn-i ^ (^aft^ ^ S^^^^o;"*-) (50) 



with 



=e«^^ =S»'' + <e (51) 

The result (50) can be easily inferred by inserting the splittings (7) and (8) into 

(49) and recalling the prescription S° = for the generators of symmetries. 
Alternatively, expression (50) can be obtained, in a clearly more involved way, 
from the variational Lagrangian 

L'iL2,E),X) = -(K/Oeaaifai...a„fa„i?"^''^ A...A^«'''"' Ai^se" (52) 

-(W0eaai6i...a„6„T<^ A R'''^' A ... A R^n-^b„.^ ^ £^^^nbn 

(obtained from the contraction of the equations (14) with the Lie derivatives of 
fields with respect to the vector field (27)) via the two-step procedure outlined 
above. 

Let us now turn the attention to the Lovelock-Chern-Simons theory de- 
scribed through the field equations = ensuing from the Lagrangian 

^According to our notation 

dx^^ A dx"! A ... A dx'''^-^ A (ix""-! A dxP = ie«flPMi -!-„-i ^g^^ (43) 
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(18). * The two step calculation of the conserved quantities starting from the 
variational Lagrangian L'(L3,^) = —^^£^gxp is obviously a rather cumber- 
some task. Nevertheless we have outlined in the introduction that conserved 
quantities can be calculated directly from (50) if wc set T = and if wo are able 
to select a suitable covariant lift of the spacetime vector field ^ up to a vector 
field S into the configuration bundle Y of the Lagrangian £2- Namely if we 
are able to select a suitable rule S"^ = $'^''(^^, rfjyC^, • . •) which determines the 
coefficients of S in terms of the components of ^ and, obviously, in terms of the 
dynamical fields; see diagram (28). As already explained the prescription which 
allows to rule out the sought for lift is the requirement that the variational La- 
grangian L'(L2,E!) and the variational Lagrangian L'(L3,^) are equivalent (so 
that conserved quantities calculated from them are equivalent, too). 

According to [3, 19, 23, 27, 28] the generalized Kosmann lift S = ^l'(^), 
obtained by the rule 

is the only one which restores the equality between the variational Lagrangians. 
This property is easily proved. Indeed, with the expression (53) we have 

•fx(€)e« = ^e^-'/e^M. (54) 



so that 



1 dL2 „ ^ 



where we have taken into account that, if T = 0, then L2{e,j^u))\j,_Q = L3{g{e)) 
and ^l^e"" = 

We also stress that the lift (53) is globally well defined. Indeed, the Kosmann 
lift -ftr(^) of a spacetime vector field ^ transforms tcnsorially (as it can be inferred 
through a direct inspection) and can be therefore globally defined. Namely, all 
local expressions (53) can be patched together to define a unique global vector 
field; see [19, 27, 28]. 

It follows that the potential for the Lovelock-Chern-Simons theory can be 
obtained directly from (50) by setting T = and taking into account the pre- 

*Wc remark that field equations are simply obtained by varying L3 with respect to the 
metric alone since the variation of the curvature terms yields a total derivative. This is the 
reason why field equations are of second order only; [39, 40]. 
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scription (53), i.e.: 

A R"-^^-^ A ... A R«-n-lbn-l 

(55) 



According to (40) the variation of conserved quantity (and therefore the vari- 
ation of energy) is then obtained from the above expression after a suitable 
integration on a spacctimc region. 

Let us remark that it is also possible to separate into (55) the contributions 
ensuing from each Lagrangian Cp of (18). Indeed, if we take the definition (15) 
into account expression (55) can be conveniently rewritten as: 

n 

U{L'{L3,i),X)=Kj2ciMp (56) 

p=0 

where 

i^P = ^^^), ^abca,b,...a„_,b„_, { (e° e"^^ V^r^D^^ + {elCnSu;'^''} 

AR"-^^^ A ... A iJ^p-i^-i A e"" A e''" A ... A e"—^ A e^"--" (57) 

and Cp arc given by (19). A direct comparison between (56) and (18) leads to 
our master formula 

U{Cp) = Up (58) 

for the covariant potential relative to the term Cp. Therefore we are now able 
to calculate the potential for a generic Lovelock Lagrangian 

n 

Llov = ^ apCp{fg) (59) 

p=0 

with completely arbitrary constant coefficients ap. It turns out to be: 

n 

Ulov = ^ apUp (60) 

In particular Cq = ^/g gives no contribution to the potential, while the contri- 
bution ensuing from the Hilbert term C,\ = ^ R is: 

i^f = V9 C^i'^'^^'^K + i^^^H") (61) 

Formula (61) reproduces exactly the formula for the variation of conserved 
quantities found in [10, 12, 26, 37] for General Relativity. This property was 
proven explicitly in [3] when the spacetime dimension is three; but the same 
calculation can be carried on, step by step, in any higher spacetime dimension. 
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Wc also point out the remarkable simple expression (57) and (60) for the po- 
tential in Lovelock gravity. In the sequel we shall test its viability, i.e. we shall 
analyse whether formula (60) allows to reproduce physically expected quantities 
in specific applications. 

Remark 2.1 Let us consider again formula (49). It is straightforward to verify 
that, on-shell, the potential satisfies: 

dU{L'{Li,E),X)=oj{Li,X,£sA) (62) 

where the {D — l)-form: 

Uj{Li,X, £sA) = nKeABCDAiBi...Ar,-iBr^-i 

is the symplectic current relative to Li; see [46]. This remarkable property 
ensures that we can make the identification 

W(L'(Li, S), X) = SxV{L,,E) - B{L,,E, X) (64) 

whore V(Li, S) is the Noether superpotential, while B{Li,E, X) is the so-called 
Reggc Teitelboim covariant correction term, namely the term which has to be 
added to the variation of the Hamiltonian to have a well posed variational 
problem (or, in other words, it is the factor which exactly cancels the boundary 
term coming from the variation of the Hamiltonian; see [22, 34, 36, 44]). This 
property, a posteriori, endows the definition we suggested for the variation of 
the energy, i.e: 

5xEj:{E)= [ U{L'{L,,E),X):^ [ u;{L,,X,£sA) (65) 
JOS is 

with a well-defined formal meaning since it is coherent with a (Hamiltonian) 
symplectic description of field theories. 

Moreover, if S is a Killing vector for the solution, i.e. £sA = 0, the sym- 
plectic form LO is vanishing, the potential U becomes a closed form and formula 
(65) is the conservation law for energy; see [3, 36, 37]. 

Inserting (7) and (8) into (62) and setting T = and S = K{^), we easily 
obtain the conservation law for Lovelock-Chern-Simons theory: 

dU{L'{L^, KiO), X) = u;{L3, X, £k(0>^, £k(0^) (66) 

where the symplectic current [47] is now given by: 

0J{L3, X, £su;, £se) = ^ ea6ca.6i...a„_.6„_i^*^^ A ... A R-n-^bn-i 

A(i:K(€)e" A Su^'' - de" A £k(0^^'') (67) 



15 



and is vanishing whenever K{^) is a KiUing vector for the dynamical fields 

Remark 2.2 The electromagnetic field can be included into the theory simply 
by minimal coupling the electromagnetic field to the Lovelock-Chern-Simons 
Lagrangian (18), i.e. 

LlZ = L3{f9) + LU9,F) (69) 
where Lem{g, F) is the Maxwell Lagrangian in (2n + 1) dimensions: 

L,m{9, F) = --J—^F^,F„0g'^"g''^ (70) 

The potential relative to the pure electromagnetic part can be calculated, through 
the two-step procedure outlined in section 2, starting from the variational La- 
grangian 

where 

(72) 

Wc shall here skip the details of the calculation since they arc just a mere and 
straightforward application of the formalism. The result turns out to be: 

^(F"^^" + 2^^"^'^^") 5Ap (73) 



f2 



D-2 



It then follows the the potential for the Lagrangian (69) is obtained from (55) 
and (73), i.e. 

U^^{L\LZ,, (€, X)),X) = U'^f{L'{Ls,0, X) + U^P{L'{L,m, (^, x)), ^) (74) 



3 The D = 5 Spherically Symmetric Charged So- 
lution 

To check the reliability of the formula (74) we test it on an exact solution and 
we show that it allows to recover the correct expression for the first law of black 

^We remark that the Lie derivative of the vielbein with respect to the Kosmann lift is given 
by (54) , while the Lie derivative of the spin connection turns out to be [23] : 

i;K(4)< = e^''e''lTi;57SM (68) 
where 7^7^ are the Christoflel symbols of the metric Q/j,^ = riabs'^si- 
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holes mechanics when stationary black holes solutions for the 5-diniensional 
Lovelock-Maxwell gravity are considered [6, 7, 11, 15, 16]. 

Let us consider a five dimensional solution of the field equations ensuing 
from the D = 5 Lagrangian (69). The gravitational contribution to the total 
gravitation-electromagnetic potential (74) can be read out from formula (56) 
with n = 2. It is given by: 

U"^{L'{Ls, 0,X)= U'^^iL'iLH, ^),X)+ U"^{L'{Lgb,0, X) (75) 

where 

W^^iL'iLH, 0, X) = ^ V5 S^J^iW-^e'^Je; + ^e^^elScof) (76) 

and: 

U^^{L'{Lgb,0,X) = ^e'^^P'^'^Cabcde Rf^ieli^S^l^ + ele'-^W^i^Sel) (77) 

are the contributions ensuing, respectively, from the Hilbert Lagrangian Ci and 
the Gauss-Bonnet Lagrangian £2. The electromagnetic contribution to the 
total potential (74) is instead given by formula (73): 

U'^P{L'{L,mAi,x)),X) = -±^5{^F^P){x + A,e) 

- ^ {F°"^e + 2^^''F^^P) SAp (78) 

ZTT 

We remark that the area O3 of the three dimensional unit sphere is 127r^ so that, 
from (5), we infer that the coupling constant k turns out to be k = 1/12-k'^. 

The spherically symmetric charged black hole solution is; see [6]: 



f\r) = l + ^^-^M + l-^ (79) 

^ A = ipdt, <P = ^ 

where d^s is the standard metric on the three-dimensional unit sphere. The 
event horizon is located at r = r+, where r+ is the positive solution of /^(r+) = 
0; see [6]. If we now consider the vector field (^,x) = {dt,Xo)t with xo = 
constant, and we integrate expression (74) on a three-sphere S of constant r 
and t (with r > r+) we obtain 

6xQs{S,0 = SM^^ (80) 
from the gravitational part (75); and 

SxQemiS, xo) = ^ + XoSQ (81) 
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from the electromagnetic part (78). The variation of the total charge is therefore 
obtained by their sum: 

SxQT::,{S,^,X)=SM + xoSQ (82) 

and reproduces the espectcd physical values. We remark that even though (80) 
and (81) depend explicitly on the radius of the sphere, their sum (82) does not 
depend on it. This property follows from the consideration we made in Remark 
2.1. Indeed from (72) we can easily infer that the vector field (^, xo) is a Killing 
vector field for the solution (79) and consequently the potential (74) is a closed 
three-form. Therefore: 



SxQToviS, Xo) = <5xQi™„(S", Xo) 



(83) 



whenever S and S' are homologous surfaces. 

In particular, we can identify Xo with (minus) the electric potential on the 
horizon xo = —f {'''+) = {^^^ thermodynamical parameter conjugated 

to the charge). In this case formula (82) becomes the first law of black hole 
thermodynamics: 

SM - ip{r+)dQ = T5S 



where: 



is the temperature T 



T=l- 

4tt 



2r. 



J_ d f\r) 
A-K dr 



of the black hole, while 



(84) 
(85) 



r=r-|- 



877 



3/2 



(86) 



is its entropy; see [6]. In obtaining this result, in the right hand side of (84) we 
have taken into account that the mass and the charge of the solution are related 
to the radius of the horizon r+ via the relation P{r+) = 0. 



4 The Entropy in Lovelock Gravity 

As it is manifest from (86) the entropy for the spherically symmetric charged 
solution differs from the area law of Einstein gravity. This is a well-known re- 
sult. It was indeed demonstrated in [35] that in Lovelock gravity each higher 
derivative term gives a contribution to the entropy in its own. We shall now 
analyse to what extent formulae (57) and (60) are in agreement with the results 
already existing in literature. 

Let us consider the Lovelock Lagrangian (59) and let us consider a stationary 
black hole solution of its field equations, admitting a Killing horizon Ti, and a 
bifurcation surface H; see [33, 35, 48, 50]. Let us denote by ^ the unique Killing 
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vector field which is null on the horizon and vanishes on H, i.e. the vector field 
which describes the evolution of a system of ROTORS (co-rotating observers; 
see [49]). We remark that on the bifurcation surface it holds true that: 

V^C, = -~2KHP^,u (87) 

where kr is the surface gravity and /S^y denotes the binormal to H. 

Let us then consider a spacelike hypersurface E which extends from spatial 
infinity to tlic^ horizon and intersects the horizon in the bifurcation surface. Since 
^ is transverse to D the variation of energy is identified with the the integral 
JgU{£_, X) of the potential W on a (2n — 1) -dimensional surface S enclosing the 
singularity. But the potential is a closed form since ^ is a Killing vector (see 
Remark 2.1). Therefore the integral JgU{^,X) is equal to Jg,U{^,X) for any 
surface S' homologous to S. In particular we can identify S' with the bifurcation 
surface H and calculate the variation of the energy on it: 

Sm = I )^llds^& = / b^fds^i, (88) 

where the potentials U^^ are given by (57). As the vector ^ vanishes on H only 
terms containing derivatives of ^ survive into (88), i.e. 

;/a/3| _ P a/3pAiii'i.../i„_ii'„_i 

V Jr"!''! DOp-ibp-i „ap bp a„_i fe„_i 1 

^ Y'-Hivi ■ ■ ■ -"-/ip-ifp-i '^IJ.p'^Vp ■ ■ ■ <^/i„_i<^L/„_i J 

Taking into account the property (87) together with the integration rule 

/ \^f9 Vl^'^ldsc^ = / v^d^-^x V["«/?„;3 (90) 
J S J s 

(which holds true for any {D — 2)-form V = \^/g V'^'^ldSa/S and any [D — 2)- 

rcgion S with induced metric g and binormal f3af3 = n^aU^j, where tIq. and U/s 
are the outward pointing unit normals to S) we obtain: 

/ l;Ufds^0 = '^l''" ^ [ V^d^-^X (91) 

7^2 f {2n-2py.2P-^ Jh^ ^ ' 

tAai/3i...ap_i/3p_iAip!^p...j/„_i-"-;ji,yi ' • Vp-i'-'p-i c"'-'p 

where gMi-;*2„-i = g,ii/^i.../i2„-i^^^ jg (2n - l)-dimensional Levi-Civita 
totally skew-symmetric tensor density (from now on a "tilde" over a quantity 
will always denote the projection of the quantity itself onto the surface H). 
Notice that in the above formula the Riemann tensors are fully projected onto 
the surface H. Moreover, since is a bifurcation surface, its extrinsic curvature 
vanishes and (via the Gauss-Codazzi equations [41, 50]) we can identify the 
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projected Riemann tensor with the {2n — 1) -dimensional Riemann tensor R 
built out of the induced (2n — l)-metric g^„, see [35]. Contracting the product 
of the Levi-Civita tensor densities in (91) and rewriting them as a product of 
(2n — 1) -dimensional Kronecker deltas we finally obtain: 

SE,{0 := T5S, = ^J^V~9 d^'-'x KZ^X-X-Mt. ' " " Kr-lt-l 

(92) 

where the temperature T is given by T = As recognized in [35] the above 
expression can be formally integrated since: 

Extracting the temperature from (92) we thence obtain: 

Sp=4np [ C(p-i){g) (94) 
Jh 

namely, the contribution to the entropy ensuing from the term Cp{g) of the 

Lovelock Lagrangian is obtained from the Lagrangian >C(p_i)(,g) defined with 
respect to the induced (2n — l)-metric 3^1/. Therefore (88) corresponds to the 
first law of black holes thermodynamics^": 

5E{£,) = T6S (95) 

where 

n „ 

S = 4n^pap / £(p_i)(5) (96) 
p=i -"^ 

This result agrees exactly with the result obtained by Wald within a Lagrangian 
approach [33] or by Jacobson- Myers within a Hamiltonian method [35]. 

We stress again that (88) is the integral of a closed form (since ^ is a Killing 
vector) and therefore the integral on H equals the integral on any other surface 
S homologous to H itself. Nevertheless on any other surface other than H the 
formal calculation of 6E{S^) = JgU{£,,X) becomes rather more involved since 
also the terms linear in ^ have to be retained and the simplifications related to 
the particular geometry of the bifurcation surface do not hold true any longer. 
Despite of this, even though the formal calculation above requires the presence 
of the bifurcation surface, on the contrary, as we did already remark in [20, 21], 
its existence become a unessential requirement in specific applications. In the 
charged spherical symmetric solution we have dealt with in section 3, e.g., the 
bifurcation surface is out of the domain of coordinates (and its presence would 
be revealed only in the max;imal analj^tic extension with Kruskal-Szekeres-like 

We remark that, according to our choice of the vector ^, the variation SE of energy turns 
out to be SE = SM — f2(') J(j) where: M is the mass of the black hole, J(i), with i = 1 . . . n, 
are the angular momenta associated with the commuting Killing vectors generating spatial 
rotations and Q^'' are the horizon angular velocities; see [16, 35]. 
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coordinates; [31, 33, 43, 48]). However the first law (84) of black hole me- 
chanics has been established performing calculation on any surface enclosing 
the singularity. This noteworthy property, we stress again, follows from the 
cohomological properties the potential inherits in presence of a Killing vector. 

5 Acknowledgments 

Wc arc grateful to A. Borowicc of the University of Wroclaw for having di- 
rected our attention to the Chern-Simons formulation of Lovelock gravity. We 
also thank L. Fatibene and M. Ferraris of the University of Torino for useful 
discussions and suggestions on the subject. 

References 

[1] A.Achucarro, P.K.Townsend, Phys.Lett. B 180 (1986), 89. 

[2] G. AUcmandi, M. Francaviglia, M. Raiteri: Class. Quant. Grav. 19 (2002), 
2633-2655 (gr-qc/0110104). 

[3] G. AUcmandi, M. Francavigha, M. Raiteri: Class. Quant. Grav. 20 (2003), 

483-506 (gr-qc/02 11098). 

[4] S.C.Anco, R.S.Timg, J. Math. Phys 43 (2002), 3984; 
S.C.Anco, R.S.Tung, J. Math. Phys 43 (2002), 5531. 

[5] M. Bahados, C. Teitelboim, J. Zanelh: Phys. Rev. Lett. 69 (1992), 1849- 
1851 (hep-th/9204099); 

M. Banados, M. Hcnncaux, C. Teitelboim, J. Zanelli: Phys. Rev. D 48 

(1993), 1506-1525 (gr-qc/9302012). 

[6] M. Bafiados, C. Teitelboim, J. Zanelh: Phys. Rev. D 49 (2) (1994), 975-986 
(gr-qc/9307033). 

[7] M. Baiiados: Phys. Rev. D 65 (2002), (hep-th/0109031); 

M. Banados: Nucl. Phys. Proc. Suppl. 88 (2000), 17-26 (hep-th/9911150). 

[8] G. Barnich, F. Brandt, Nucl.Phys. B633 (2002) 3 (hep-th/0111246). 

[9] N.D. Birrel, P. C. W. Davies: Quantum fields in curved space, Cambridge 
University Press (1982). 

[10] I. Booth: PhD Thesis from University of Waterloo (gr-qc/0008030); 
I. Booth, R.B. Mann: Phys. Rev. D 59, 064021 (gr-qc/9810009); 
I. Booth, R.B. Mann: Phys. Rev. D 124009 (1999) (gr-qc/9907072). 

[11] D.G. Boulware, S. Deser: Phys. Rev. Lett. 55 (24) (1985), 2656-2660. 



21 



[12] J. D. Brown, J. W. York, Phys. Rev. D 47 (4), (1993), 1407; 
J. D. Brown, J. W. York, Phys. Rev. D 47 (4) (1993), 1420. 

[13] A. H. Chamscddinc: Nucl. Phys. B 346 (1990), 213-234; 

A. H. Chamscddinc, J. Frohhch: Comm. Math. Phys. 147 (1992), 549. 

[14] C. M. Chen, J. M. Nester: Gravitation & Cosmology 6, 257, (2000) (gr- 
qc/0001088). 

[15] J. Crisostomo, R. Troncoso, J. ZaneUi: Phys.Rev. D 62 (2000), 084013 
(hep-th/0003271). 

[16] M. H. Dehghani: Phys. Rev. D 67 (2003), 064017 (hep-th/0211191). 
[17] S. Deser, B. Tekin: Phys. Rev. D 67 (2003), 084009 (hep-th/0212292). 

[18] L. Fatibene and M. PrancavigHa, Natural and gauge natural formalism for 

classical field theories: a geometric perspective including spinors and gauge 
fields, in Pubbhcation for Kluwcr Academic Publishers. 

[19] L. Fatibene, M. Ferraris, M. Francaviglia, M. Godina: Gen. Rel. Grav. 30 

(1998) , 1371-1389. 

[20] L. Fatibene, M. Ferraris, M. Francaviglia, M. Raiteri: Annals of Phys. 275 

(1999) , 27. 

[21] L. Fatibene, M. Ferraris, M. Francaviglia, M. Raiteri, Phys. Rev. D60, 

124012 (1999); 

L. Fatibene, M. Ferraris, M. Francaviglia, M. Raiteri, Phys. Rev. D60, 

124013 (1999); 

L. Fatibene, M. Ferraris, M. FrancavigUa, M. Raiteri, Annals of Phys. 284, 
197 (2000), gr-qc/9906114. 

[22] M. Ferraris and M. FrancavigUa: Atti Scm. Mat. Univ. Modena, 37 (1989), 
61; M. Ferraris and M. Francaviglia: Gen. Rel. Grav., 22, (9) (1990), 
965; M. Ferraris, M. Francaviglia, I. Sinicco: // Nuovo Cimento, 107B, N. 
11, (1992), 1303; M. Ferraris, M. Francaviglia: 7th Italian Conference on 
General Relativity and Gravitational Physics, Rapallo (Genoa), September 
3-6, 1986. 

[23] M. Ferraris, M. Francaviglia, M. Raiteri: Class. Quant. Grav. in press (gr- 
qc/0305047). 

[24] M. Ferraris, M. Francaviglia, in: Mechanics, Analysis and Geometry: 200 
Years after Lagrange, Editor: M. Francaviglia, Elsevier Science Publishers 

B. V., (Amsterdam, 1991) 451; 

A. Trautman, in: Gravitation: An Introduction to Current Research, L. 
Witten ed. (Wiley, New York, 1962) 168; 

A. Trautman, Commun. Math. Phys., 6, (1967), 248. 



22 



[25] M. Ferraris, M. Francaviglia and 0. Robutti, in: Geometrie et Physique, 
Proceedings of the Journees Relativistes 1985 (Marseille, 1985), 112 - 125; 
Y. Choquet-Bruhat, B. Coll, R. Kerner, A. Lichnerowicz eds. (Hermann, 
Paris, 1987). 

[26] M. Francaviglia, M. Raiteri: Class. Quant. Grav. 19, (2002), 237-258 (gr- 

qc/0107074). 

[27] M. Godina, P. Mattcucci: in press for J. Geom. Phys. (matli.DG/0201235); 

P. Mattcucci, gr-qc/0201G79. 

[28] M. Godina, P. Matteucci, J. A. Vickers: J. Geom. Phys. 39 (2001), 265-275. 

[29] H.Goldschmidt, D.Spencer, J. Diff. Geom. 13 (1978) 455. 

[30] M.B. Greens, J.H. Schwarz, E.Witten: Superstring Theory (Cambridge 
University Press, England 1987); 

D. Lust, S. Theusen: Lectures on String Theory (Springer, Berlin, 1989); 
J.Polchinski: String Theory (Cambridge University Press, England 1998). 

[31] S. W. Hawking, G.F.R. Elhs, The Large Scale Structure of Space-Time, 
(Cambridge University Press, Cambridge 1973). 

[32] S. W. Hawking, C. J. Hunter, Class. Quant. Grav. 13, (1996) 2735; 
S. W. Hawking, G. T. Horowitz, Class. Quant. Grav. 13, (1996) 1487. 

[33] V. Iyer and R. Wald: Phys. Rev. D 50 (1994), 846 (gr-qc/9403028). 

[34] V. Iyer, R. M. Wald: Phys.Rev. D 52 (1995), 4430-4439 (gr-qc/9503052). 

[35] T.Jacobson, R. Myers: Phys. Rev. Lett. 70, 3684-3687 (1993) (hep- 
th/9305016). 

[36] B. Julia, S. Silva: Class. Quant. Grav. 15, (1998), 2173 - 2215 (gr- 

qc/9804029); 

B. Julia, S. Silva: Class. Quant. Grav. 17 (2000), 4733 - 4744 (gr- 
qc/0005127). 

[37] J. Kijowski: Gen. Rel. Grav. 29 (1997), 307. 

[38] I. Kolaf, P.W. Michor, J. Slovak: Natural Operations in Differential Ge- 
ometry, Springer- Verlag (New York, 1993). 

[39] C. Lanczos: Ann. Math. 39 (1938), 842. 

[40] D. Lovelock: J. Math. Phys. 12 (1971), 498; 
D. Lovelock: J. Math. Phys. 13 (1972), 874. 

[41] C.W. Misner, K.S. Thorne, J.A. Wheeler, Gravitation (Freeman, San Fran- 
cisco, 1973). 



23 



[42] R. C. Myers: Phys. Rev. D 36 (1987), 392. 

[43] I. Racz, R. M. Wald: (gr-qc/9507055). 

[44] T. Regge, C. Teitelboim: Annals of Physics 88 (1974), 286. 

[45] G. Sardanashvily: Generalized Hamiltonian Formalism for Field Theory, 
World Scientific (Singapore, 1995). 

[46] S. Silva, Nucl. Phys. B 558, (1999), 391 (hep-th/9809109). 

[47] C. Teitelboim, J. Zanelli: in Constraints Theory and Relativistic Dynam- 
ics, Proceedings of the Workshop held in Florence, Arcetri, Italy (World 
Scientific, Singapore 1987), 94. 

[48] P. K. Townscnd: Black Holes, DAMPT University of Cambridge. 

[49] M. Visscr: Phys. Rev. D 48 (1993), 5697-5705 (hep-th/9307194). 

[50] R. M. Wald: General Relativity, University of Chicago Press (Chicago, 
1984). 

[51] E. Witten: Nucl. Phys. B 311 (1988), 96; 
E. Witten: Nucl. Phys. B 323 (1989), 113. 



24 



